Concatenating Decoherence Free Subspaces with Quantum Error Correcting Codes 
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An operator sum representation is derived for a decoherence- free subspace (DFS) and used to 
(i) show that DFSs are the class of quantum error correcting codes (QECCs) with fixed, unitary 
recovery operators, and (ii) find explicit representations for the Kraus operators of collective deco- 
herence. We demonstrate how this can be used to construct a concatenated DFS-QECC code which 
protects against collective decoherence perturbed by independent decoherence. The code yields an 
error threshold which depends only on the perturbing independent decoherence rate. 
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Decoherence-free subspaces (DFSs) have recently 
emerged as an alternative way to protect fragile 

quantum states against decoherence, alongside "conven- 
tional" quantum error correcting codes (QECCs) |@,|| 
and the new "dynamical decoupling" schemes p]| . This 
is of particular importance in quantum computation, 
where the promise of a speed-up compared to classical 
computers hinges crucially on the possibility to maintain 
quantum coherence throughout the computation fl^| . So 
far, DFSs and QECCs have been considered as distinct 
methods, often characterized as "passive" and "active" 
respectively. However, as we will show here, in fact DFSs 
can be considered as a special class of QECCs, character- 
ized as having a particularly simple form of recovery op- 
erators. Conditions for the existence of non-trivial DFSs 
are stringent: the decoherence process should be "col- 
lective", meaning that the bath couples in a symmetric 
way to all qubits. So far conditions for collective de- 
coherence have been formulated in a Hamiltonian form 
HQ , and in the Lindblad semigroup form ||U . Here we 
will present an alternative formulation in terms of the 
operator sum representation (OSR) jDj], which has the 
advantage of establishing a direct link to the theory of 
QECCs. This OSR formulation enables us to combine 
DFSs and QECCs, into a concatenated scheme which 
can error-correct the more general physical situation of 
"cluster decoherence" . For the price of longer codewords, 
this concatenated scheme operates with a substantially 
reduced error threshold. 

Hamiltonian Formulation of Decoherence Free 
Subspaces. — Conditions for DFSs within the general 
(non-Markovian) framework of semigroup dynamics were 
derived in Ref. M. We first briefly rederive these con- 
ditions in a simplified form. Consider a closed quantum 
system, composed of a system S of interest defined on 
a Hilbert space Ti (e.g., a quantum computer) and a 
bath B. The evolution of the closed system is given 
by psB{t) — Up5s(0)U^, where the unitary evolution 
operator (we set h = 1) is U = cxp(— iHt). The full 
Hamiltonian is H = Hg <8)Ib+Is®Hb + Hj, where 
Hs, Hg and H/ are, respectively, the system, bath and 
interaction Hamiltonians, and I is the identity opera- 
tor. Assuming initial decoupling between system and 



bath, the evolution of the closed system is given by: 
PSB(t) = U[ ( os(0)<X>/9b(0)]U' 1 '. Quite generally, the inter- 
action Hamiltonian can be written as H/ = £ Q F Q <g> B Q , 
where F Q and B Q are, respectively, system and bath op- 
erators. Suppose that there exists a degenerate subset 
{|fc)} of eigenvectors of the F Q 's such that: 



F a \k)=a a \k) Va,\k). 



(1) 



If Hs leaves the Hilbert subspace Ti =Sp[{\k)}] invari- 
ant, and if we start within Ti, then the evolution of the 
system will be decoherence free (DF). To show this, ex- 
pand the initial density matrices of the system and the 
bath in their respective bases: ps(0) = £ 
/>b(0) =£^^|/i)(z/|. Using Eq. 
combined operation of the bath and interaction Hamil- 
tonians over Ti as: 
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one can write the 
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This clearly commutes with Hs over Ti. Thus since 
neither Hs (by our own stipulation) nor the combined 
Hamiltonian H c takes states out of the subspace: 



U[|?> ® |m>] 



u 



'U c |Ai>, 



(2) 



where XJx = cxp(— iHxt), X = S,c. Hence it is clear, 
given the initially decoupled state of the density ma- 
trix, that the evolution of the closed system will be: 



PSB{t) = EijSyU 



>S\i 
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follows using simple algebra that after tracing over the 
bath: p s (t) = Tr B [p SB {t)} = U S( os(0)U s , i.e., that the 
system evolves in a completely unitary fashion on Ti.: un- 
der the condition of Eq. (Q) the subspace is DF. As shown 
in Ref. U, Eq. H is also a necessary condition for a DFS. 

Operator Sum Representation on a Decoherence Free 
Subspace. — In the OSR, the evolution of the density 
matrix is written as: ps(t) =Trs[U(ps ® /°s)U^] = 
£ a A a ps(0) AJ, where the "Kraus operators" are given 
by: 



A a = vM/'lUl 



(3) 



(|/x), \v) are bath states) and satisfy the normalization 
constraint J2 a A l A a = Is 
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Let Ti. be an TV-dimensional DFS. In this case it follows 
immediately from Eqs. (§) and (|J) that the Kraus oper- 
ators all have the following representation (in the basis 
where the first N states span H): 



9aUs 
A a 



g a = Mv\VcW). (4) 



Here A a is an arbitrary matrix that acts on "H (H = 
H(B'H ± ) and may cause decoherence there; Us is Us re- 
stricted to TL. This simple condition can be summarized 
as follows: 

Theorem I. A subspace Ti is a DFS iff all Kraus op- 
erators have an identical unitary representation upon re- 
striction to it, up to a multiplicative constant. 

Thus, in the OSR, the task of identifying a DFS re- 
duces to finding a subspace in which all the Kraus oper- 
ators act as the system unitary evolution operator. We 
now give an example for the important case of collec- 
tive decoherence (CD). CD is generally described within 
the following scenario: the system operators {F Q } in 
the interaction Hamiltonian form the Lie algebra su(2) 
HH]. This means that the interaction Hamiltonian can 
be rewritten as: 



V, 



S_ ® V_ 



(5) 



Here S Q = X^iLi °~? are global Pauli spin operators (i is 
the qubit index) satisfying the sl(2) commutation rela- 
tions, and V a are the bath operators coupled to these 
degrees of freedom. A more restricted case of CD arises 
when only phase damping processes are allowed, so that 
V + = V = 0. We will concentrate on this case, as it is 
fully analytically solvable. For simplicity we will assume 
throughout that Hg = 0. 

Pure phase damping on a single qubit is described by 
the Pauli <r z matrix. In the standard basis it is easy 
to verify that the matrix representation of the global 
phase damping operator S z is: S z = diag[/(j)], where 
f(j) = (no. of 0's) - (no. of l's) in the binary rep- 
resentation of j (j = 0...2 A_1 ). E.g., for two qubits: 
S z = diag[2, 0, 0, —2]. Since S z is diagonal, the action of 
the interaction Hamiltonian H/ = S z (£> V z can be writ- 
ten simply as: H/|j)|f) = f(j)\j) <8> V^f). Hence the 



action of the full Hamiltonian is: H|j) 



\J) 



$7, 



where = f(j)V z + H B . Similarly, = \j) ® 

whence exp(— iHt)\j) \v) = \j) ® exp(— iftjt)\v), 
so that the Kraus operators can be evaluated explic- 
itly: (j'\A a \j) — 6jf^/v(ii\ exp(— iftjt)\v). Thus, the 
Kraus operators for pure phase damping CD have a 
diagonal matrix representation in the standard basis, 
which can be written compactly as: A a = dia; 

V^(mI ex P(~ i$ljt)\v). For example, in the case 
of two qubits we obtain A a = diag 



o (2) o (0) 
ya , ga ■ 



o (0) 
ya , 



,(-2) 



with ga = exp(— iHgi)]^) and 

y/v(n\ exp(-i[H B ± 2\ z ]t)\v). By Theorem I the three 



(±2) 
9 a 



blocks in A a correspond to three DFSs. Next, consider 
the effect of this OSR on a general density matrix (omit- 
ting standard algebra): [p(0)] jk i — > J2a [ A ap(°) A I] ifc = 

\PWtikL,»v "9^ ■ Qiiv ■ A s expected, no mixing of 
density matrix elements occurs. The time depen- 
dence of each element is determined by the sum in 
the last expression, which, motivated by the under- 
standing that decoherence is taking place, we write 
formally as a decaying exponential (although without 
a Markovian approximation Poincare recurrences may 

occur). Thus: exp(-t/r jk ) = Y,y, v v 9ifj 9$ = 
£^|exp{^[/(j)-/(fc)](V. + [V Z ,H B ]) + 
with higher order terms, all depending on powers of 
f(j) — f(k), given by the Campbell-Hausdorff formula. 
In agreement with the general theory of DFSs , the 
DF states are those for which f(j) = f(k), in which case 
l/ r ife = 0- Normally (i.e., in the "pointer basis" [||), all 
other states are expected to have 1/Vjfc > 0, although in 
order to verify this one must specify V z . This expecta- 
tion is confirmed for a harmonic bath |Lj . 

Decoherence Free Subspaces as Quantum Error Cor- 
recting Codes. — Quantum error correction can be re- 
garded as the theory of reversal of quantum operations on 
a subspace |l0| |. This subspace, C = Sp[{|ii)}], is inter- 
preted as a "code" (with codewords {|«l)}) which can be 
used to protect part of the system Hilbert space against 
decoherence (or "errors") caused by the interaction be- 
tween system and bath. The errors are represented by 
the Kraus operators {A a } ||. To decode the quantum 
information after the action of the bath, one introduces 
"recovery" operators {R r }- A QECC is a subspace C 
and a set of recovery operators {R r }. Ref. || gives two 
equivalent criteria for the general condition for QECC. 
It is possible to correct the errors induced by a given set 
of Kraus operators {A a }, (i) iff 



■RVA a — 











Vr, a, 



or equivalently, (ii) iff 



(6) 



(7) 



In both conditions the first block acts on C; B,. a and A a 
are arbitrary matrices acting on C 1 (Ti = C ® C ). Let 
us now explore the relation between DFSs and QECCs. 
First of all, it is immediate that DFSs are indeed a valid 
QECC. For, given the (DFS-) representation of A Q as in 
Eq. (|j), it follows that Eq. (]?]) is satisfied with "f ao = g*gt- 
Note, however, that unlike the general QECC case which 
has a full-rank matrix 7 a &, in the DFS case this matrix 
has rank I (since the a th row equals row I upon mul- 
tiplication by gi/gT), implying that a DFS is a highly 
degenerate QECC 

A DFS is an unusual QECC in another way: decoher- 
ence does not affect a perfect DFS at all. Since they are 
based on a perturbative treatment, other QECCs (e.g., 
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stabilizer or GF(4) codes @) are specifically constructed 
to improve the fidelity to a given order in the error rate, 
which therefore always allows for some residual decoher- 
ence to take place. The absence of decoherence to any 
order for a perfect DFS is due to the existence of sym- 
metries in the system-bath coupling which allow for an 
exact treatment. These symmetries are ignored by per- 
turbative QECCs either for the sake of generality, or be- 
cause they simply do not exist, as in the case of inde- 
pendent couplings. Given a DFS, the only "errors" that 
can take place involve the unitary rotations of codewords 
(basis states {\i}} of Ti) inside the DFS, due to the sys- 
tem Hamiltonian Hg (this may actually be the desired 
evolution if one is implementing a computation inside Ti 
using Hs). Thus, the complete characterization of DFSs 
as a QECC is given by the following: 

Theorem II. Let C be a QECC for error operators 
{A a }, with recovery operators {R r }. Then C is a DFS 
iff upon restriction to C, R r oc Ijjy for all r. 

Proof. First suppose C is a DFS. Then by Eqs. (|J) and 
(|), 



( 9atJ S 


1 




\ o 







To satisfy this equation, it must be true that 

R = f^S C r \ 

The condition g a UjjD r = implies D r = by unitar- 
ity of Us. Also, since A a is arbitrary, generically the 
condition C r A a = implies C r = 0. Thus upon re- 
striction to C = Ti, indeed R r oc U s (by unitarity of 
Us, \Ka/g a \ = !)• Now suppose R r oc TJ' S . The very 
same argument applied to A a in Eq. (JsJ) yields A a oc Us 
upon restriction to C. Since this is exactly the condition 
defining a DFS in Eq. (||), the theorem is proved. 

We conclude that DFSs are a particularly simple in- 
stance of general QECCs, where upon restriction to the 
code subspace, all recovery operators are proportional to 
the inverse of the system evolution operator. 

Quantum Error Correction On A Decoherence Free 
Subspace. — A non-ideal DFS will still be subject to some 
decoherence Q. DFSs are efficient under conditions in 
which each qubit couples to the same environment (col- 
lective decoherence). Ordinary QECCs are designed to 
be efficient when each individual qubit couples to a dif- 
ferent environment (independent decoherence). While 
neither code is efficient in the extreme limit when the 
other is, QECCs will still work for correlated errors 
whereas DFSs will not work in the independent error 
case Q . One would generally expect the likelihood of K 
qubits collectively coupling to the same environment to 
decrease with increasing K 0. Thus an interesting sit- 
uation ("cluster decoherence") arises when small blocks 
of qubits undergo collective decoherence (e.g., groups of 



neighboring identical atoms on a polymer chain), while 
this symmetry is broken perturbatively by independent 
decoherence between blocks. Here we show how by 
adding an additional layer of QECC encoding, the DFS 
can be stabilized against such computational errors. 

In the collective decoherence case [Eq. (0)], the small- 
est DFS which can encode one logical qubit is made up of 
four physical qubits || . Consider an operator basis which 
covers all possible errors which can occur on this 4-qubit 
DFS. From the discussion above, we know that a DFS has 
Kraus operators which all are direct sums of a fixed uni- 
tary transformation on Ti and variable transformations 
(whose exact form is irrelevant) on Ti . Thus, pertur- 
bative errors (of size e) on a DFS can conveniently be 
represented by Kraus operators with the following struc- 
ture: 

Here A a represents the (dominant) contribution due to 
ideal collective decoherence, and the second matrix rep- 
resents the symmetry breaking perturbation. Qi (a 2 
x 2 matrix for the 4-qubit case) acts just on the DFS; 
Q2 (2 x 14) takes states from Ti into Ti, but we need 
not worry about this as a separate process, since it can 
also be corrected using QECC inside Ti; Q3 (14 x 2) 
takes states from Ti into Ti ± ; Q4 acts just on Ti and 
is irrelevant to our discussion. Thus, to first order in 
time, all of the relevant errors can be enumerated as: (i) 
Independent errors acting on the encoded DFS states. 
A basis for these errors are the Pauli operators: X, Z, 
Y = XZ, and I acting only on the DFS qubits |0l) and 
(e.g., X|0i) = (ii) Errors which take the 

system into Ti 1 - = Sp[{|ji)}^£ 2 ]. Define 14 operators 

Pj = \jl)({0l\ + (U|), Pj : H h 7Y-L, out of a total of 
28 in the Q3 block. In order to cover all possible errors 
which might take the DF states outside of Ti, it suffices to 
consider the effect of Pj and PjZ, where Z is the phase 
error operator acting on Ti as defined above. The possible 
errors which can occur on our 4-qubit DFS arc thus given 
by: £ = {X, Y, Z, P 3 , PjZ}. The task is now to find an 
appropriate error correction scheme. To do so, we may 
use the DF states |0l) and | ljr ) to construct the well- 
known "perfect" 5-qubit QECC juj. This concatenation 
yields new encoded states |0b) and composed of 20 
physical qubits. In the standard approach to correcting 
errors with the 5 qubit code, one uses a quantum net- 
work to calculate the syndrome on ancilla qubits, and 
uses this syndrome to apply the appropriate correction 
procedure on the encoded qubits. Thus, in order to apply 
a similar procedure to correct for the standard bit flip er- 
rors which now act on our 4-qubit DF states, we simply 
convert the gates used in the standard error correction 
procedure to gates which perform the same operations 
on |0i) and |1l) and which do not disturb the states in 
Ti . Under these conditions, it is obvious that X, Y, Z 
can be corrected if the DF qubits decohere independently 
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(the condition on standard QECC), as we stipulated. 
Furthermore, consider the following modified controlled- 
not gate: C|0 L ,0 L ) = |0 L ,0 L ), C|1 L ,0 L ) = |1 L ,0 L ), 
C|j'l,0l) = plus unspecified operations (which 

do not have to concern us) which ensure that C is uni- 
tary. Using this gate, one can correct for the effect of Pj 
on the system as follows: attach an ancilla DF qubit in 
the \0l) state, perform C on this state, and repeat the 
procedure on every DF qubit of the 5-qubit codeword. 
This detects whether one of the Pj errors has occurred. 
Essential to this procedure is the fact that we do not 
disturb the original |0l) and states. Once an er- 
ror of this type has been detected, one may recover the 
erred state to |0l). The error has then been reduced 
to a standard Pauli one on a DFS state, and can be 
fixed by QECC. For example (the first qubit, x — or 1, 
belongs to the codeword and the second is the ancilla): 

\xl)\0l) S {\x L ) + \2 L ))\0 L ) £ \x L )\0 L ) + \2 L )\2 L ). Next 
the ancilla is measured: if the result is 0, no error has oc- 
cured; if it is 2 then the state is recovered to \0l) and 
the resulting standard Pauli error can be fixed by com- 
parison of the erred DF qubit to the other 4 DF qubits 
in the codeword. Finally, the extreme case of single in- 
dependent physical qubit errors is also dealt with by the 
concatenated DFS-QECC code, because the above proce- 
dure automatically also corrects the errors which are the 
subject of the standard QECC procedures, representing 
them as linear combinations drawn from our natural DFS 
basis of errors £. Of course, the 20-qubit code presented 
above is less efficient than the standard 5-qubit QECC 
code if these independent errors are the dominant ones. 

The final question regarding the concatenated DFS- 
QECC scheme concerns the threshold for fault tolerant 
quantum computation |Te| ] . The threshold probability of 
error has been estimated to be on the order of 10 -6 per 
operation In the present context it is simplest to 

discuss this issue in the language of the Markovian semi- 
group master equation. Following the notation of Ref. || , 
we consider error generators {F Q } yielding a decohering 
term [F a , p(t)¥t] in the master equation. Assume these 
generators produce errors with a rate A. Then the fidelity 
F(t) — Tr[ps(0)ps(t)] is generally reduced by a term of 
0(X 2 t). Upon inclusion of a symmetry-breaking pertur- 
bation by error generators {G p } of order e, with A» t, 
we find, following the arguments in Ref. ||, that this 
leads to a decrease in the fidelity, not of order O(Xet) as 
one might naively expect (since [G p , p(t)Ft] terms have 
appeared), but, remarkably, only of order 0(e 2 t) (due to 
the [G p ,p(t)G\} terms) ||. 

Perturbing a DFS thus produces error rates which are 
solely proportional to the strength of the additional per- 
turbing process, and which are not dependent on the rate 



of the error process which generates the DFS. Applying 
a QECC to this system, as discussed above, improves 
the fidelity to 1 — 0(e A t 2 ). Now in quantum computa- 
tion one generally envisions a realization of a quantum 
computer which has an a priori low decoherence rate A. 
But in the ordinary QECC case, it is this rate A which 
sets the error threshold. Hence, in the scenario envi- 
sioned here, we have effectively decoupled the rate A from 
the error threshold, and a significant improvement may 
therefore result, if the additional perturbation e is suffi- 
ciently small. Concatenation of a DFS with QECC thus 
has the potential to achieve the goal of truly fault toler- 
ant quantum computation, not just quantum memory. 
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